Symmetry, Integrability and Geometry: Methods and Applications SIGMA 3 (2007), 065, 11 pages 

The Rahman Polynomials Are BispectralB 

F. Alberto GRUNBAUM 

Department of Mathematics, University of California, Berkeley, CA 94720, USA 



E-mail: grunbaum@math. berkeley. edu 



Received February 01, 2007, in final form April 22, 2007; Published online May 03, 2007 



Original article is available at http://www.emis.de/journals/SIGMA/2007/065/ 



i Abstract. In a very recent paper, M. Rahman introduced a remarkable family of polyno- 

■ mials in two variables as the eigenfunctions of the transition matrix for a nontrivial Markov 

Cn| I chain due to M. Hoare and M. Rahman. I indicate here that these polynomials are bispectral. 

This should be just one of the many remarkable properties enjoyed by these polynomials. 
d \ For several challenges, including finding a general proof of some of the facts displayed here 

the reader should look at the last section of this paper. 

ff) ' Key words: bispectral property; multivariablc polynomials; rings of commuting difference 

operators 

<^ ■ 2000 Mathematics Subject Classification: 33C45; 22E45 

' Dedicated to the memory of Vadim Kutznetsov 

a : 

I first met Vadim at a meeting in Esterel, Canada back in 1994. I only saw him a few more 
times after that, in Amsterdam, in Leeds, in Louvain-la-Neuve, and the last time in Leganes. 

y—i • Although we never wrote any papers together we shared a number of common interests. I always 

found Vadim extremely helpful and willing to share his ideas and his vast knowledge. He was 
' also willing to put up with people whose work was not in the mainstream and he usually would 

have some useful remark or a suggestion to make. I like to think that he would have enjoyed 
seeing this paper which builds on some remarkable work of another common friend I met for the 

JO | first time at the same Esterel meeting. 

o : 

^ ! 1 Preliminaries 

> ■ 

Starting about thirty years ago, M. Hoare and M. Rahman, published a number of papers 
studying a class of statistical models that were nicely associated with "classical" orthogonal 
■ - - ■ polynomials in one variable. The reader can consult [T9j l20l [2Tj EE]. 

In a very recent paper [22], which I have seen before publication thanks to the kindness 
of M. Rahman, they take a huge step and consider the case of several variables. I propose 
to call these statistical models with the name Hoare-Rahman, and to call the multivariable 
polynomials that arise with the name Rahman polynomials. I am thankful to Prof. Rahman for 
several inspiring conversations on the contents of [22], in different coffee houses and restaurants 
in Fremont, California on the last days of 2006. 

The model of interest here captures very well the idea of a game of chance where the player 
takes a risk that can never improve his/her winnings but is then given a second opportunity 
to improve his/her lot. I am convinced that models of this kind should be of great interest 
in several areas of applied mathematics not only in physics (where they were conceived) but 
also in biology, population dynamics, evolutionary models, etc. On top of their potential for 
applications the polynomials discovered by M. Rahman are a piece of beauty. I hope that this 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 



2 



F.A. Griinbaum 



paper will incite other people to look around for the gems contained in [22]. I only recall the 
points from this paper that are needed in the present context, urging the reader to go back to 
the source for inspiration. 

For a systematic study of polynomials in several variables one should consult [8] as well as 
the work on Macdonald polynomials of various kinds, see [31]. A look at the pioneering work of 
Tom Koornwinder, see for instance [28], is always a very good idea. 

2 The Rahman polynomials 

Given generic values of the parameters px, P2, Ps, Pi (these could even be complex numbers) 
define t, u, v, w by means of 

{Pi + P2){Pl + Pz) „. _ (Pi + Pz){PA +P3) 



Pi {pi +P2+P3+ Pa) ' Pz{pi +P2+PZ+ Pi) ' 

(Pi +P2){P2 + Pa) [PA +P2){PA +P-3.) 

w - 



P2(pi +P2+P3 + Pa) ' Pa(pi +P2+P3+ Pa) ' 

Now we can define with M. Rahman the polynomials in x, y given by 



P ( \ - srsrsrsr (- m )i+j(- n )k+e(-x) i+k {-y)j+e 

P mAx , V)-^^ i\j!kW(-N) i+j+k+e 



% v?v k w l . 



I will not try to explain M. Rahman's path to these polynomials, but the reader should con- 
sult [22]. In particular one finds in [22] a discussion of the connection with the 9 — j symbols 
used in quantum mechanics. 

3 The Hoare— Rahman model and the kernel /C 

Denote by 

b(k,N,a) = r^Ja k (l-a) N - k 

the usual binomial distribution giving the probability of k successes in N repetitions of an ex- 
periment with two possible outcomes: success (with probability a) and failure (with probability 
1 — a). Likewise denote by 62(^1,^2! N; 0i, P2) the trinomial distribution with parameters j3\, 
This means that that for i±, 12 nonnegative and such that i\ + 12 < iV we have 



b 2 (ii,i2,N; fix, fa) = ^ a (l -fh- P 2 ) N ~ ll - t2 



ii\i 2 \{N — %x — %i)\ ) ' 



Define now a transition probability kernel 

K-(jx,j2;k,i2) 

in the space of non-negative integer physical variables i\, 12 and ji, 22 such that 

h + i 2 < N, 31 +J2<N 
by means of the expression 

min(iijl) min(j 2 J2) 

^(jl, h'lh.h) = ^2 XT b(ki,i 1 ;a 1 )b(k2,i2;a 2 ) 

fc 1= fc 2 =0 
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x h(ji - h,32 - k 2 ,N - ki - k 2 ;Pi,P 2 )- 



The expression above is written so as to agree with the one in [22]. Notice that its interpre- 
tation is that of the probability of a transition from the state (ii,i 2 ) to the new sate (ji,j2) in 
one step (consisting of three tosses) of our process, as will be explained in the next section. 

At this point the parameters an, a 2 , Pi, P2 are free to lie anywhere in the interval [0, 1], with 
the sum of Pi and 02 not exceeding the value 1. 

The corresponding Markov chain is not a "nearest neighbour one" in the sense that the 
classical Ehrenfest urn model, see \IU[ [231 132j . or all other classical models are. Here the 
transition probability matrix has all its entries positive and therefore there is no sense in which 
a one-step transition moves you only to your immediate neighbours. I find it remarkable that 
the complete spectral problem can be solved for such a kernel, but this is exactly what the 
Rahman polynomials do, at least under certain extra conditions. 

4 Poker dice 

A concrete example of the model is provided by a gambling game known as Poker dice. The 
version considered by M. Hoare and M. Rahman goes as follows. You have N identical dice 
with n faces each. There are three kinds of faces: red, black and "blanks". The probabilities of 
getting a red or a black face when rolling one dice are a\ and a 2 , respectively. The state of the 
system is given by (11,12), the number of red and black faces among the iV dice. Think of the 
remaining N — i\ — i 2 dice as failures. 

Start form the state (11,12) and get ready to roll three times. The resulting number of reds 
and blacks after these three rools will be denoted by (j'1,,7'2). This will constitute one time step 
in the Markov chain to be described now. 

The ii dice that show red are tossed and yield ki red faces. The i 2 dice that showed a black 
face are tossed and yield fc 2 blacks. Now you take the initial N — %i — i 2 blanks and combine 
them with the ii — ki dice that failed to reproduce their original red color and the i 2 — k 2 dice 
that failed to reproduce their black color. The total number of "failures" is now N — ki — k 2 . 
These dice are thrown by themselves: they yield a certain number of red and black faces which 
are added to the ki red faced dice and fc 2 black faced dice that were saved earlier. This total 
number of reds is now ji and the total number of blacks is j 2 . To obtain full generality the 
probabilities of getting a red or a black for each individual dice on this last toss can be taken to 
be Pi and p 2 , not necessarily the old values ai and a 2 . 

The case of iV = 5 is referred to as standard "poker dice" in [22] and will be used to illustrate 
later on the main point in this paper, namely the bispectral property of the Rahman polynomials. 
We have verified this property for many other values of N, but we do not have a nice general 
expression for the operator in question, or a nice proof of its existence for general N. 

5 The orthogonality weight 

In [22] one sees that the Rahman polynomials are orthogonal with respect to the trinomial 
distribution 

b 2 (x,y;N, 971, 172), 

where the parameters rji, 772 are given in terms of pi, P2, P3, Pa by means of 



P1P2O1 +P2 + P3 + P\) 



V2 = 



P3Pi(Pl +P2+P3+P4) 



m = 



(Pi +P2)(Pl +P3)(P2+Pi) 



(Pi + P3)(P4 + P2)(P4 + P3) ' 
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Given this weight, for an arbitrary choice of the parameters rji, r] 2 , one could blindly decide 
to "order" the set of monomials 

x l y j 

and then apply the Gram-Schmidt process to obtain orthogonal polynomials. Of course there 
is no unique way to order these monomials and each ordering will produce a different family of 
polynomials. 

One of the features of the work of M. Rahman is that it automatically gives one set of 
polynomials that are already orthogonal and, as we see in [22] and here, enjoy remarkable 
properties. Moreover he produces a four parameter family of these polynomials, having started 
from a the two parameters 771, ij 2 . There is clearly a lot of inspiration behind this construction, 
and the reader should consult [22j . 



6 The bispectral property 

Recall that KL depends on one set of four parameters and the Rahman polynomials on a different 
set of four parameters. 

As shown in |22j . the Rahman polynomials satisfy the eigenvalue relation 



for (x,y) in "physical space", i.e. x, y are integers satisfying x > 0, y > 0, x + y < N, with 
j 1 - air/i - a 2 r?2 - r?it(l - Qi) - r/ 2 u(l - cx 2 ) 



1 — airji — a 2 r] 2 
1 - air/i - a 2 r/2 - r?ii>(l - «i) - r] 2 w(l - a 2 ) 1 " 
1 - a\r}x - a 2 r/2 / 
= (1 - - (3 2 u) m (l - 0iv - fow) n . 

More explicitly, with the definitions made above this reads, 

ji ,32 

For this to be true we need to impose some relations connecting the parameters pi, p 2 , P3, P4 
and a±, a 2 , 0i, (3 2 . This issue has been addressed in [22], 

In the considerations above we are thinking of these polynomials as functions of the physical 
space variables x, y for fixed values of the frequency parameters m, n. 

In the spirit of [5] one could ask if they are the eigenvectors of another spectral problem 
when considered, for fixed x, y, as functions of m, n. If this were true we would have an analog 
of the three-term recursion relation satisfied by orthogonal polynomials of one variable. In the 
continuous-continuous case this question was raised in [5] for a second order Schrodinger type 
differential operator. For many ramifications of this problem the reader can see |18| . 

We give below the simplest recursion that we have been able to find. To be very concrete we 
display the result in the case of = 5 which corresponds to the Poker dice example described 
above. The operator in question consists of multiplication by the function (p±+p 2 )x— {p^+PA)y 
and it gives a "five-point" formula, where most points in the "spectral space" have four neigh- 
bours. There are, as indicated above, similar recursions for other values of N but we do not 
have a general expression that would cover all cases. 

In the case of iV = 5 the frequency space (and also the physical space) looks as follows 
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If we declare two points to be adjacent when they can be obtained by going one step either 
in the vertical or horizontal directions, and we order the 21 states as follows 

(0,0), (1,0), (2,0), (3, 0), (4, 0), (5, 0), (0, 1), (1, 1), (2, 1), (3, 1), (4, 1), 
(0, 2), (1, 2), (2, 2), (3, 2), (0, 3), (1, 3), (2, 3), (0, 4), (1,4), (0, 5) 

then one can work out easily the adjacency matrix. 

The matrix B to be detailed below will have nonzero entries exactly where the adjacency 
matrix has a 1. It is therefore as much of a local matrix as one can get in this setup. It will 
correspond to the analog of the three term recursion relation of the usual situation when one 
deals with orthogonal polynomials in one variable. 

We claim that the Rahman polynomials, as functions of m, n satisfy the difference equation 

BP m ,n(x,y) = ((pi +P2)X - (P3+Pi)y)Pm,n( X ^y) 

as long as (x, y) lie in physical space. 

Here the matrix B is made up of rows as specified one at a time below. Notice that in each 
row there is an element denoted by Si for which we do not give an expression. This is not needed 
since the sum of all the entries in each row equals the value zero. 

[si, a u 0, 0, 0, 0, h , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

with 

= 5fflP 3 (P4 +P2)(P4 +P3+P2+Pl) b = 5p 2 (P3 + Pl)P4{P4 + P3 + P2 + Pi) , 

(P3 + Pi ) (PlPi ~ P2P3 ) (P4 + P2 ) {PlPA ~ P2P3 ) 

[a 2 , 52, 02, 0, 0, 0, 0, da, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

with 

_ PlP4 ~ P2P3 _ 4piP3 (g4 + P2 ) (P4 + P3 + P2 + Pi ) 

P3+P1 ' {P3 + Pl){PlPA ~ P2P3 ) 

d _ 4p 2 (P3 +Pl)Pi(P4: +P3 +P2 
(P4 +P2){PlP4 -P2P3) 

[0, a 3 ,s 3 , c 3 , 0, 0, 0, 0, da, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

with 

_ 2(pip 4 -P2P3) _ 3piff 3 (P4 + P2)(P4 + P3 + P2 + Pi) 

P3+P1 (P3+Pl) (PlPA ~ P2P3 ) 
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d _ 3ff 2 (j>3 + Pl)P4(P4 + P3 + P2 + Pi) , 
{PA+P2)(PlPA-P2P3) 

[0, 0, a 4 , s 4 , c 4 , 0, 0, 0, 0, d 4 , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 



with 



3(pip 4 - P2P3) 2pip 3 (p 4 + p 2 ){p4 +P3+P2+P1) 

a 4 = , C4 - 



d 



P3+PI (P3 + Pi ) (PlP4 - P2f»3 ) 

2p2(P3 + Pl)P4(P4 + P3 + P2 + Pi) 



{Pl+P2)(piPi -P2P3) 

[0,0,0,a 5 ,S5,C5,0,0,0,0,e 5 ,0,0,0,0,0,0,0,0,0,0] 



with 



4 (piP4 - P2P3 ) PlP3 (P4 + P2 ) (P4 + P3 + P2 + Pi ) 

05 — ; , c 5 



e 5 



P3 + Pi (P3 + Pl)(PlP4 - P2P3) 

P2(P3 +Pl)P4(P4 +P3 +P2 +Pl) 



(P4 +P2)(PlP4 -P2P3) 

[0, 0, 0, 0, a 6 , s 6 , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

with 

5(PiP4 -P2P3) 

«6 = ; ; 

P3 +Pl 

[dr, 0, 0, 0, 0, 0, s 7 ,c 7 , 0, 0, 0, a 7 , 0, 0, 0, 0, 0, 0, 0, 0, 0] 

with 

4j>2(P3 + Pi)Pi{P4 +P3+P2+P1) 
(P4 +P2XP1P4 -P2P3) 

4 +P2)(P4 +P3 +P2 +. 
3 +Pl)(PlP4 -P2P3) 

[0, a 8 , 0, 0, 0, 0, 6 8 , s 8 , d 8 , 0, 0, 0, e 8 , 0, 0, 0, 0, 0, 0, 0, 0] 



_ 4pip 3 (p 4 + £» 2 )(P4 + P3 + P2 + Pi) ^ _ PlP4 ~ P2P3 , 

(P3 + Pi) (P1P4 - P2P3 ) P4 + P2 



with 



PlP4 - P2P3 , PlP4 - P2P3 

as = ; , k 



d 



Pa +P2 P3+ Pi 

3piP3(P4 + P2)(P4 + P3 + P2 + Pi) 3p 2 (P3 + Pl)P4(P4 + P3 + P2 + Pi) 



(P3 + Pl)(PlP4 - P2P3) (P4 + P2)(P1P4 - P2P3) 

[0, 0, e 9 , 0, 0, 0, 0, c 9 , sg, a 9 , 0, 0, 0, d 9 , 0, 0, 0, 0, 0, 0, 0] 



with 



_ 2pi£> 3 Q 4 + P2 ) (P4 + P3 + P2 + Pi ) _ 2Qip 4 -P2P3) 
(P3 + Pi ) (P1P4 - P2P3 ) P3 + Pi 

d = 2P2 (P3 + Pi )P4 (P4 + P3 + P2 + Pi ) _ PlP4 ~ P2P3 , 

(P4 + P2) (PlP4 - P2P3) P4 + P2 

[0, 0, 0, dm, 0, 0, 0, 0, 010, 810, cio, 0, 0, 0, do, 0, 0, 0, 0, 0, 0] 



with 



_ 3(pip 4 -P2P3) _ PlP3(P4 +P2)(P4 +P3 +P2 +Pl) 

P3 + Pi (P3 + Pi ) (P1P4 - P2P3 ) 
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d _ PlP4 - P2P3 e _ P2 (g3 + Pi )P4 (g4 +P3+P2+Pl) , 

P4+P2 ' (PA+P2)(PlP4~ P2P3 ) 

[0, 0, 0, 0, 6n , 0, 0, 0, 0, a n , an, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 



with 



4(piP4 - P2P3) , PiP4 ~ P2P3 
an — ; > oil — ; ! 

P3 +Pl P4+ P2 

[0, 0, 0, 0, 0, 0, ai 2 , 0, 0, 0, 0, s 12 , d 12 , 0, 0, c 12 , 0, 0, 0, 0, 0] 



with 



2(piP4 - 3p 2 (P3 + Pl)P4(P4 +P3+P2+Pl) 

ai2 = ; , C12 - 



d 



12 



£>4 + P2 (P4 + P2){PlP4 ~ P2P3) 

3P1P3(P4 + P2)(P4 +P3+P2+Pl) 



{P3+Pl)(PlP4 -P2P3) 

[0, 0, 0, 0, 0, 0, 0, C13 , 0, 0, 0, 013, sis, di 3 , 0, 0, e 13 , 0, 0, 0, 0] 



with 



PlP4 - P2P3 2(PlP4 - P2P3) 

a i3 — ; , C13 — ■ , 

P3 +Pl P4+ P2 

2pip 3 (p 4 + p 2 )(p4 +P3+P2+ Pi) 2p 2 (p3 + Pl)P4{P4 + P3 + P2 + Pi) 

Ct]_3 — ^13 — ' 

(P3 + Pi ) (PlP4 - P2P3 ) (P4 + £»2 ) (P1P4 - P2P3 ) 

[0, 0, 0, 0, 0, 0, 0, 0, d u , 0, 0, 0, oi 4j si4, cm, 0, 0, e u ,0, 0, 0] 



with 



_ 2(pip 4 ~ P2P3) _ PlP3 {P4 + P2) (P4 +P3+P2+Pl ) 

P3+P1 (P3+Pl) (PlP4 ~ P2P3 ) 

d _ 2(PlP4 -P2P3) _ P2 (P3 + Pi )P4 (P4 +P3+P2+Pl) . 

P4 + P2 (P4 + P2 ) (PlP4 ~ P2P3 ) 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 015 , 0, 0, 0, C15 , sis, 0, 0, 0, 0, 0, 0] 



with 



2 {piP4 ~ P2P3 ) 3 (pip 4 - p 2 p 3 ) 
«15 — ; , C15 — ■ ; 

P4 +P2 P3+ Pi 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, oie , 0, 0, 0, sis, cis, 0, die, 0, 0] 



with 



3(j>lP4 -P2P3) _ 1piP3ij>4 + P2)ij>4 + P3 + P2 + Pi) 

P4+P2 {P3+Pl) (PlP4 ~ P2P3 ) 

d _ 2p 2 (P3 +Pl)p4(P4 +P3 +P2 + Pl) . 
16 {P4 +P2){PlP4 -P2P3) 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ei 7 , 0, 0, a 17 , s m c 17 , 0, d 17 , 0] 



with 



_ PlP4 ~ P2P3 _ PlP3 (g4 + P2) (g4 + P3 + P2 + Pi ) 

P3 + Pi (P3 + Pi ) (PlP4 - P2P3 ) 

d _ P2(P3 +Pl)P4(P4+P3 +P2+Pl) _ 3(pip 4 -P2P3) , 

(P4 + P2 ) (P1P4 - P2P3 ) P4 + P2 
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[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ois , 0, 0, cis, ais, 0, 0, 0] 

with 

3 (PIPA ~ P2P3 ) 2 (ptpi - p 2 P3 ) 

Ol8 — ; , Ci8 ~ 



Pi +P2 P3+ Pi 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, c 19 , 0, 0, s 19 , 019, d 19 ] 

with 

PlP3(P4 +P2)(P4 +P3 +P2+Pl) 



«19 



{P3 +Pl)(piP4 -P2P3) 



4(PlP4 - P2P3) , P2{P3 + Pi)Pa{PA +P3+P2+P1) 

C19 = ; j "19 



Pi + P2 ijPA + P2 ) (P1P4 - ) 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, a 20 , 0, c 20 , s 20 , 0] 



with 



4(pi£>4 - P2P3) P\Pa ~ P2P3 

°20 — ; , C20 — ; ; 

Pa +P2 P3+ Pi 

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, a 2 i , 0, s 21 ] 

with 

5(piP4 -P2P3) 



"21 



PA+P2 



We note that these relations hold for arbitrary (generic) values of pi, p 2 , P3, p±. 

If one considers multiplication by x or by y we get a "seven-point" formula. Of course it 
would be interesting, in the spirit of [3] to consider the full algebra of difference operators going 
along with the Rahman polynomials. 

Thinking along these lines one could also consider the full algebra of difference operators that 
go along with the Rahman polynomials when thought of as functions of x, y. This simply means 
the algebra of all matrices that commute with the kernel KL. While the original kernel K, has 
a very nice probabilistic interpretation there is no reason why other operators in this algebra 
should allow such an interpretation. More on this in the next section. 



7 Global vs. local 

One should naturally think of the kernel K. given earlier as a global operator, one that in this 
discrete case is represented by a full matrix and that in the continuous case would be represented 
by an integral operator. At the core of many considerations leading up to the bispectral problem 
is the issue of the existence of a local operator that would commute with the global one. This 
remarkable piece of work is forever associated with the names of D. Slepian, H. Landau and 
H. Pollack. This is discussed for instance in |13|. [14] . 

Here, as an illustration we show that this question has a positive answer in the setup of Hoare 
and Rahman. Take N = 3 in which case we can see, for arbitrary choice of the free parameters 
a±, a 2 , /?i, fa that the matrix given below is essentially the only one that commutes with the 
kernel fC. Once again, we have checked that this holds for other values of N, and this case is 
used here only for illustration purposes. 

In the spirit of the work just alluded to one should consider /C as the analog of the operator 
of time-band-time limiting and the operator below as the analog of the second order differential 
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operator that comes from separation of variables for the Laplacian in R 3 in prolate spheroidal 
coordinates. 

In agreement with the case N = 5 considered as an illustration earlier, we order the states 
in physical space as follows 

(0, 0), (1, 0), (2, 0), (3, 0), (0, 1), (1, 1), (2, 1), (0, 2), (1, 2), (0, 3). 

Then, up to trivial scaling and shifts the only local matrix that commutes with /C is given by 





X2 








x 3 














" 


X 4 


X5 


xe 








x 7 

















X8 


xg 


xw 








Xn 

















X12 


X\3 































Xl5 


X16 





xn 











a; is 








X19 


X20 


X21 





X22 











X23 








X2A 


X25 























X26 








X27 


£28 


X29 

















£30 





X3I 


^32 


























X33 





X34. 



as long as we choose its entries as follows 

xi = 3(aia 2 /?2 - 2ai/3 2 + aia 2 ft - a 2 0\ - ol\0\ - a\a 2 + ot\)/ (ctift), 



X2 



3a 2 /9i 



X3 



X4 



(ai-l)a 2 (/? 2 +01-1) 



aift aift 
x 5 = {2a\a 2 2 + ot 2 2 ~ 5u 1 2 + 2a 1 a 2 /3 1 - a 2 0i - 3ai0 1 - 2a\a 2 - a 2 + 3ai)/ (aift) 
2a 2 ft _ 2(ai - l)a 2 (ft + ft - 1) 



x 7 



x 8 



aift ' aift 
29 = (aia 2 p 2 + 2a 2 /3 2 - 4aift + «ia 2 ft + a 2 ft - 3aift - aia 2 - 2a 2 + 3ai)/(aift), 
a 2 ft 1 3(«i - l)a 2 (ft + ft — 1) 



aJio 



Z13 



3>16 



^19 



^21 
^24 

x 2 e 



X 2 g 



Xu = 1, Xl 2 
aift 

3(a 2 -ai)(ft + ft - 1) 



xu 



aift 

(a 2 -l)(ft + ft-l) 



aift ft 
X15 = 2(aia 2 /3 2 - 2aift + ai« 2 ft - a 2 ft - aift - aia 2 + ai)/(aift), 



xn 



2a 2 ft 
aift ' 
(ai - l)a 2 (ft + ft - 1) 
aift 



£18 



(a 2 -l)(/?2 + fti-l) 
ft 



£20 = (aia 2 ft + a 2 ft - 3aift + aia 2 ft - 2aift - aia 2 - a 2 + 2ai)/ (aift 



a 2 ft 



£22 = 1, 



aift 

2(ai - l)a 2 (ft +ft - 1 
aift 

2(a 2 -l)(ft+ft-l) 



£23 



(a 2 -l)(ft+ft -1) 



£25 



ft 

2(a 2 - ai )(/3 2 + ft-l) 
aift 



x 2 7 = (aia 2 (3 2 - 2aift + a\a 2 0\ - a 2 ft - aift - ct\a 2 + ai)/(aift), 



Q201 
Oil 02 



X29 



h 



X30 



2(a 2 -l)(ft + ft-l) 

ft 
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X31 = 



(ai - l)a 2 (/3 2 + /3i - 1) 
3(a 2 -l)(/3 2 +/?i-l) 



a^3 2 = 



(a 2 -«i)(/? 2 + 
ai/3 2 



^33 = 



X34 = 0. 



8 Some closing challenges 

Here is a challenge that is motivated by recent work on matrix valued orthogonal polynomials. 
This subject was initiated by M.G. Krein a long time ago, see [29l [30]. There has been quite 
a bit of work in terms of producing explicit examples of these polynomials. For a representative 
sample see [3j [6j [Til CG3 EH EH ESI El] • Most of this work deals with looking for examples 
that would enjoy the appropriate bispectral property. On the surface of it, there is little relation 
between these polynomials and the topics discussed here. But appearances can be deceiving. 
In a recent paper [12] I noticed how the example in [11] has a probabilistic interpretation (still 
waiting to be unravelled) and leads to nonnearest neighbours transitions on the nonnegative 
integers. This goes beyond the usual birth-and-death processes for which there is a well known 
relation to orthogonal polynomials, see [U [H [231 12H ESI [27]. In the same paper [12] one finds 
examples arising in the theory of matrix valued orthogonal polynomials (in one variable) that 
give rise to Markov chains in two dimensional discrete space. 

Another challenge would be to find the appropriate concrete algebraic geometry that lies 
behind the relation among the variables a±, a 2 , (3±, /? 2 and pi, p 2 , ps, p^. In particular this 
could have some bearing on the situation discussed in [22] where the spectrum of fC is degenerate, 
see (5.15). 

One challenge mentioned already a couple of times consists in writing down the general 
formula for a five-term recursion formula that is displayed here only for N = 5. 
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